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& Abstract 

We study a class of nonstandard interactions of the newly discovered 125 GeV Higgs- like res- 
onance that are especially interesting probes of new physics: flavor violating Higgs couplings to 
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leptons and quarks. These interaction can arise in many frameworks of new physics at the elec- 

D 

troweak scale such as two Higgs doublet models, extra dimensions, or models of compositeness. We 
derive constraints on flavor violating Higgs couplings using data on rare decays, electric and mag- 

> 

netic dipole moments, and meson oscillations. We find that flavor violating Higgs boson decays to 

ON 

leptons can be sizeable with, e.g., h — >■ tu and h — > re branching ratios of 0(10%) perfectly allowed 
by low energy constraints. We estimate the current LHC limits on h — > Tfi and h — > re decays 

o 

by recasting existing searches for the SM Higgs in the rr channel and find that these bounds are 
already stronger than those from rare tau decays. We also show that these limits can be improved 
^ significantly with dedicated searches and we outline a possible search strategy. Flavor violating 

Higgs decays therefore present an opportunity for discovery of new physics which in some cases 
may be easier to access experimentally than flavor conserving deviations from the Standard Model 
Higgs framework. 
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I. INTRODUCTION 



Both ATLAS and CMS have recently announced the discovery of a Higgs-like resonance 
with a mass of m h ~ 125 GeV [1-4], further supported by combined Tevatron data [5]. 
An interesting question is whether the properties of this resonance are consistent with the 
Standard Model (SM) Higgs boson. Deviations from the SM predictions could point to the 
existence of a secondary mechanism of electroweak symmetry breaking or to other types 
of new physics not too far above the electroweak scale. While there is a large ongoing 
experimental effort to measure precisely the decay rates into the channels that dominate 
for the SM Higgs, it is equally important to search for Higgs decays into channels that are 
subdominant or absent in the SM. For instance, since the couplings of the Higgs boson to 
quarks of the first two generations and to leptons are suppressed by small Yukawa couplings 
in the SM, new physics contributions can easily dominate over the SM predictions. Another 
possibility, and the main topic of this paper, is flavor violating (FV) Higgs decays, for 
instance into r/i or fie final states. The study of FV couplings of the Higgs boson has a 
long history [6-22], however, many of the results in our paper are new. We discuss in detail 
possible search strategies for FV Higgs decays at the LHC and derive for the first time limits 
from LHC data. 

As we will show below the constraints on many FV Higgs decays are rather weak. In fact, 
for h — » r/i and h — >■ re already now 1 , without targeted searches, the LHC is placing limits 
that are comparable to or even stronger than those from rare r decays. As we shall see later, 
re-casting a h — > tt analysis with 4.7 ftr 1 of 7 TeV ATLAS data [23] gives a bound on the 
branching fraction of the Higgs into t/j, or re around 10%. We will also demonstrate that 
dedicated searches can be much more sensitive. 

These decays could thus give a striking signature of new physics at the LHC, and we 
strongly encourage our experimental colleagues to include them in their searches. Another 
experimentally interesting set of decay channels are flavor conserving decays to the first two 
generations, e.g., h — > on which we will comment further below. We emphasize that 

large deviations from the SM do not require very exotic flavor structures. A branching ratio 
for h — ?■ Tfi comparable to the one for h — » tt, or a h — > ji + fj,~ branching ratio a few 

1 By h — > Tfi we always mean the sum of h — > t + and h — > t~ fi + and similarly for the other decay 
modes. 
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times larger than in the SM can arise in many models of flavor (for instance in models with 
continuous and/or discrete flavor symmetries [24], or in Randall-Sundrum models [25]) as 
long as there is new physics at the electroweak scale and not just the SM. The lepton flavor 
violating decay h — > t/j, has been studied in [11], and it was found that the branching ratio 
for this decay can be up to 10% in certain Two Higgs Doublet Models (2HDMs). 

In fact, there may already be experimental hints that the Higgs couplings to fermions 
may not be SM-like. For instance, the BaBar collaboration recently announced a 3.4<j 
indication of flavor universality violation in b — > ctu transitions [26] , which can be explained 
for instance by an extended Higgs sector with nontrivial flavor structure [27]. 

The paper is organized as follows. In Sec. II we introduce the theoretical framework we 
will use to parameterize the flavor violating decays of the Higgs. In Sec. Ill we derive bounds 
on flavor violating Higgs couplings to leptons and translate these bounds into limits on the 
Higgs decay branching fractions to the various flavor violating final states. In Sec. IV we 
do the same for flavor violating couplings to quarks. We shall see that decays of the Higgs 
to T/i and to re with sizeable branching fractions are allowed, and that also flavor violating 
couplings of the Higgs to top quarks are only weakly constrained. Motivated by this we 
turn to the LHC in Section V and estimate the current bounds on Higgs decays to r/i and 
re using data from an existing h — > rr search. We also discuss a strategy for a dedicated 
h — >■ Tfi search and comment on differences with the SM h — > tt searches. We will see 
that the LHC can make significant further progress in probing the Higgs' flavor violating 
parameters space with existing data. We conclude in Section VI. In the appendices, we give 
more details on the calculation of constraints from low-energy observables. 

II. THE FRAMEWORK 

After electroweak symmetry breaking (EWSB) the fermionic mass terms and the cou- 
plings of the Higgs boson to fermion pairs in the mass basis are in general 

c Y = -mJir R - y, J (.rji 1 )h + h.c. + • • • , (i) 

where ellipses denote nonrenormalizable couplings involving more than one Higgs field oper- 
ator. In our notation, = qL, are SU (2)l doublets, fn = ur, cIr, vr, £r the weak singlets, 
and indices run over generations and fermion flavors (quarks and leptons) with summation 
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implicitly understood. In the SM the Higgs couplings are diagonal, Yy = (m,i/v)Sij, but 
in general NP models the structure of the Yy can be very different. Note that we use the 
normalization v = 246 GeV here. The goal of the paper is to set bounds on and identify 
interesting channels for Higgs decays at the LHC. Throughout we will assume that the Higgs 
is the only additional degree of freedom with mass (9(100 GeV) and that the Y^'s are the 
only source of flavor violation. These assumptions are not necessarily valid in general, but 
will be a good approximation in many important classes of new physics frameworks. Let 
us now show how Yj? ^ (mi/v)Sij can arise in two qualitatively different categories of NP 
models. 

a. A single Higgs theory. Let us first explore the possibility that the Higgs is the only 
field that causes EWSB. For simplicity let us also assume that at energies below ~ 200 GeV 
the spectrum consists solely of the SM particles: three generations of quarks and leptons, 
the SM gauge bosons and the Higgs at 125 GeV. Additional heavy fields (e.g. scalar or 
fermionic partners which address the hierarchy problem) can be integrated out, so that we 
can work in effective field theory (EFT) — the effective Standard Model. In addition to the 
SM Lagrangian 

Csm = RilPfi + f R 0f R ~ [\jUlf R )H + h.c] + D^D^H - \ H (tfH - , (2) 

there are then also higher dimensional terms due to the heavy degrees of freedom that were 
integrated out: 

AC Y = -^(fif R )H(H^H) + h.c. + • • • , (3) 

Here we have written out explicitly only the terms that modify the Yukawa interactions. 
We can truncate the expansion after the terms of dimension 6, since these already suffice to 
completely decouple the values of the fermion masses from the values of fermion couplings 
to the Higgs boson (see also [15]). Additional dimension 6 operators involving derivatives 
include 

A£d = §(flYfL)(HW,H) + §(f R Yf R )(HW,H) + -.., (4) 

where (HHl%H) = HUD^H - (iD^W)H. The couplings are complex in general, 
while the \ l [ R are real. The derivative couplings do not give rise to fermion-fermion-Higgs 
couplings after EWSB and are irrelevant for our analysis. In Eq. (4) there are in principle 
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also terms of the form (f l LR ilpf LR )H^H, which, however, can be shown to be equivalent to 
(3) by using equations of motion. 

After electroweak symmetry breaking (EWSB) and diagonalization of the mass matrices, 
one obtains the Yukawa Lagrangian in Eq. (1), with 



V2m = V L 



v 2 ' 



V R v, V2Y = V L 



v 2 



Vl (5) 



where the unitary matrices Vl, Vr are those which diagonalize the mass matrix, and v = 
246 GeV. In the mass basis we can write 

where A = VlX'Vr. In the limit A — > oo one obtains the SM, where the Yukawa matrix Y is 
diagonal, Yv = m. For A of the order of the electroweak scale, on the other hand, the mass 
matrix and the couplings of the Higgs to fermions can be very different as A is in principle 
an arbitrary non-diagonal matrix. 

Taking the off diagonal Yukawa couplings nonzero can come with a theoretical price. 
Consider, for instance, a two flavor mass matrix involving r and ji. If the off-diagonal entries 
are very large the mass spectrum is generically not hierarchical. A hierarchical spectrum 
would require a delicate cancellation among the various terms in Eq. (5). Tuning is avoided 
if [28] 

with similar conditions for the other off diagonal elements. Even though we will keep this 
condition in the back of our minds, we will not restrict the parameter space to fulfill it. 

b. Models with several sources of EWSB: Let us now discuss the case where the Higgs 
at 125 GeV is not the only scalar that breaks electroweak symmetry. The modification of 
the above discussion is straightforward. The additional sources of EWSB are assumed to 
be heavy and can thus still be integrated out. Their EWSB effects can be described by a 
spurion x that formally transforms under electroweak global symmetry and then obtains 
a vacuum expectation value (vev), which breaks the electroweak symmetry. If x has the 
quantum numbers (2,1/2) under SU{2) L x U(l)y it can contribute to quark and lepton 
masses. 2 This allows the Yukawa interactions Y of the 125 GeV Higgs to be misaligned with 



A spurion which transforms as a triplet can also contribute to Majorana masses for neutrinos. 
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respect to the fermion mass matrix m in Eq. (1). 

The simplest example for a full theory of this class is a type III two Higgs doublet model 
(2HDM) where both Higgses obtain a vev and couple to fermions. In the full theory both 
of the scalars then have a Lagrangian of the form (1) 



where the index a runs over all the scalars (with imaginary for pseudoscalars), and mi 
receives contributions from both vevs. In addition there is also a scalar potential which 
mixes the two Higgses. Diagonalizing the Higgs mass matrix then also changes Yfi, but 
removes the Higgs mixing. For our purposes it is simplest to work in the Higgs mass basis. 
All the results for a single Higgs are then trivially modified, replacing our final expressions 
below by a sum over several Higgses. For a large mass gap, where only one Higgs is light, the 
contributions from the heavier Higgs are power suppressed, unless its flavor violating Yukawa 
couplings are parametrically larger than those of the light Higgs. The contributions from 
the heavy Higgs correspond to the higher dimensional operators discussed in the previous 
paragraph. This example can be trivially generalized to models with many Higgs doublets. 

We next derive constraints on flavor violating Higgs couplings and work out the allowed 
branching fractions for flavor violation Higgs decays. In placing the bounds we will neglect 
the FV contributions of the remaining states in the full theory. Our bounds thus apply 
barring cancellations with these other terms. 

III. LEPTONIC FLAVOR VIOLATING HIGGS DECAYS 

The FV decays h — > e/i, er, //r arise at tree level from the assumed flavor violating Yukawa 
interactions, Eq. (1), where the relevant terms are explicitly 

C Y D - Y efl e L fi R h - Y^ e Ji L e R h - Y eT e L r R h - Y Te f L e R h - Y^ T Ji L T R h - Y Tfl f L fi R h + h.c. . 



The bounds on the FV Yukawa couplings are collected in Table I, where for simplicity of 
presentation the flavor diagonal muon and tau Yukawa couplings, 



Cy = - 



mfif R -Y«(fif R )h a + h.c. + --- , 



(8) 



(9) 



C Y D -Y^HLHiih - Y TT r L T R h + h.c. , 



(10) 
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Channel 


Coupling 


Bound 


\x — > e7 




< 3.6 x 10~ 6 


\x — > 3e 


VIVI 2 + I^I 2 


< 0.31 


electron g — 2 


Re(^v) 


-0.019... 0.026 


electron EDM 


|im(y eM y^ e )| 


< 9.8 x 10~ 8 


\x — > e conversion 




< 4.6 x 10~ 5 


M-M oscillations 


\y 4- y* I 


< 0.079 


r — >■ e7 


Vl^el 2 + |yr| 2 


< 0.014 


r — >■ e/x/x 


Vl^el 2 + |F er | 2 


< 0.66 


electron g — 2 


Re(y er y re ) 


[-2.1... 2.9] x 10~ 3 


electron EDM 


|im(y er y re )| 


< 1.1 x 10~ 8 


r — >■ 1U7 


Vl^i 2 + ivi 2 


< 1.6 x 10~ 2 


r — > 3fi 




< 0.52 


muon g — 2 


Re(y Mr y rM ) 


(2.7 ±0.75) x 10~ 3 


muon EDM 


im(y MT y TM ) 


-0.8... 1.0 


/U — > e7 


fly y I 2 + ly y \ 2 ) 1/4 

W 1 t/j, 1 re\ ~\ 1 fiT J -eT\ ) 


< 3.4 x 10~ 4 



Table I: Constraints on flavor violating Higgs couplings to e, fi, r for a Higgs mass nih = 125 GeV 
and assuming that the flavor diagonal Yukawa couplings equal the SM values (see text for details). 
For the muon magnetic dipole moment we show the value of the couplings required to explain the 
observed Aa M (if this is used only as an upper bound one has ^Re(y Mr y rAt ) < 0.065 at 95%CL). 

were set equal to their respective SM values (Ypn) SM = m ^/v, (^rr) SM — rn T /v. Similar 
bounds on FV Higgs couplings to quarks are collected in Table II. 

We first give more details on how the bounds in Tables I and II were obtained and then 
move on to predictions for the allowed sizes of the FV Higgs decays. 

A. Constraints from r — > ^7, r — > and fj, — > e"f 

The effective Lagrangian for the r — >■ /I'j decay is given by 

£ e ff = c L Q Ll + c R Q Rl + h.c. , (11) 
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Figure 1: Diagrams contributing to the flavor violating decay r — > /ry, mediated by a Higgs boson 

with flavor violating Yukawa couplings. 

where the dim-5 electromagnetic penguin operators are 

Ql 7 ,r 7 = ^m T (fx a a ^P LiR r)F af3 , (12) 

with a, (3 the Lorentz indices and F a/3 the electromagnetic field strength tensor. The Wilson 
coefficients cl and cr receive contributions from the two 1-loop diagrams shown in Fig. 1 
(with the first one dominant), and a comparable contribution from Barr-Zee type 2-loop 
diagrams, see Fig. 12 in Appendix A. The complete one loop and two loop expressions are 
given in Appendix A. 

In the approximation <C Y TT , only the first of the one- loop diagrams in Fig. 1 is 
relevant (in addition to the 2-loop diagrams). Using also m M <C m T <C m h and assuming 
y^fi, Ytt to be real, the expressions for the one-loop Wilson coefficients cl and cr simplify 
to 



The 2-loop contributions are numerically 

<zr = >;;(-o.o82y« + oad—^ = o.^—^, m 

where in the last step we used for the top Yukawa coupling Y tt = {Y tt )sM = fn t /v = 0.67, 
and we have normalized the results to for easier comparison. (By fh t , we denote the top 
quark mass parameter in the MS renormalization scheme, rn t ~ 164 GeV.) The analytical 
form of the Wilson coefficient can be found in Appendix A. The same result applies to 
c r°° P with the replacement Y* — > Y^ T . The 2-loop contribution with the top quark in the 
loop cancels to a good approximation with the W contribution. The end result is thus very 
sensitive to the precise value of Y tt . For Y tt ~ fh t /v the 2-loop contribution is about four 
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times as large as the 1-loop contribution, while for other values of Y tt (e.g., Y tt ~ —m t /v) 
the 2-loop contribution can be an order of magnitude larger. 

In terms of the Wilson coefficients Cl and Cr, the rate for r — > /ry is 

r <^^) = ^(N 2 + W 2 )- < 15 > 

Using a Higgs mass rrih = 125 GeV and assuming Y TT = m T /v, Y tt = m t /v, we can 
then translate the experimental bound BR(r — > /ry) < 4.4 x 10~ 8 [29] into a constraint 
\/\Y TfM \ 2 + |F MT | 2 < 1.6 x 10~ 2 (see Table I). The bound is relaxed if Y TT and/or Y tt are 
smaller than their SM values. 

The expressions for // — >■ e^y and r — > are obtained in an analogous way with the 
obvious replacements (r — > fi, /j, — > e for the first and — >■ e for the second in Eqs. (12), 
(13), (14), (15)). The experimental bound BR(pi ->■ e7) < 2.4 • 1(T 12 [29] then trans- 
lates to ^[Y^l 2 + |Y^| 2 < 3.6 x 10" 6 and BR(r ->■ ery) < 3.3 x 1(T 8 [29] to a constraint 
y/\Y Te \ 2 + |^| 2 < 1.4 x 10- 2 using the SM values for Y TT , F M/1 , and Y tt . The — >■ bound 
is completely dominated by the two loop contribution, while for r — > ej, the two loop and 
one loop contributions are comparable. 

The decay /i — > cy can also be used to place a bound on the combination Y^ T Y re using 
the 1-loop Wilson coefficient 

4 loop ~ -L ^y; t y; £ (- 3 + 2 log 4 V (16) 

L 8m 2 h ^ T Te \ m 2 T J 

As before, c^ oop is obtained by replacing Y*Y* e by Y Tfl Y eT . The 2-loop contribution is propor- 
tional to Y^e and Y efl . Setting them to zero, one obtains a bound ([l^Fi-el 2 + [Y^Fei-l 2 ) 1 ^ 4 < 
3.4 x 10" 4 . 

B. Constraints from r — > 3/x, r — > e/x/x, /x — >■ 3e 

The tree level diagram in Fig. 2, with a single insertion of the flavor violating Yukawa 
coupling Y TfM , can generate the decay r — > 3/i. The decay width is given by 

r , x 5m^|y w | 2 (ivi 2 + |y;,| 2 ) 

r(T 3/i) = 2WK ' 

where we neglected the terms additionally suppressed by the muon mass. The experimental 
bound BR(r 3/i) < 2.1 x 1(T 8 [30], translates into a constraint sJ\Y?^ + \Y^ T | 2 < 0.52 




Figure 2: Tree level Higgs exchange diagram leading to r — > 3/x decay. 



for rrih = 125 GeV and = m^/v. Even though the decay r — > 3/i is tree level it leads 
to a weaker limit on Y TfM , Y^ T than r — > /i7 because of the smallness of F MAt . The same is 
true for constraints on F Me , Y efl , Y Te , Y er following from the processes /i — > 3e and r — > e/i/i, 
which are suppressed by Y ee and Y^ respectively. The bounds in Table I are obtained using 
the experimental results BR(/i -> 3e) < 1.0 x 10~ 12 [31] and BR(r e/i^u) < 2.7 x 10~ 8 [32] 
with obvious replacements for the Yukawa couplings in Eq. (17) (for r — > efifi Eq. (17) 
needs to be also multiplied by 4/5 ). 



C. Constraints from muonium antimuonium oscillations 

A ji + e~ bound state (called muonium M) can oscillate into an e + /i~ bound state 
(antimuonium M) through the diagram in Fig. 3. The time- integrated M — >■ M con- 
version probability is constrained by the MACS experiment at PSI [33] to be below 
P(M — > M) < 8.3 x 10~ n / Sbi where the correction factor S B < 1 accounts for the splitting 
of muonium states in the magnetic field of the detector. It depends on the Lorentz structure 
of the conversion operator and varies between S B = 0.35 for (S ± P) x (S ± P) operators 
and S B = 0.9 for P x P operators [33]. Conservatively, we use the smallest value S B = 0.35 
throughout. Since we will find that M-M oscillation constraints are much weaker than 
those from from /i — > cy and ji — > e conversion, this approximation suffices for illustrative 
purposes. 

The theoretical prediction for the M — > M conversion rate is governed by the mixing 
matrix element (see, e.g., [34]) 

M -/t I t I + ^ ePfi)e] + ^ ePfi)e] I t I lt\ Mtt 

J V ^MM — \ I \i.\e — I e | ^2 I I e ^fi ~ ^elfl / , 
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Figure 4: A diagram contributing to the anomalous magnetic moment g — 2 of the muon through 
FV couplings of the Higgs to r/x. 

where '[x and ix are the spin orientations of particle X. We can work in the non- 
relativistic limit here. For a contact interaction, the spatial wave function of muonium, 
0i s = exp(— r/a M )/[7ra| f ] 1 / 2 only needs to be evaluated at the origin. (Here r is the 
electron-antimuon distance and om = (m e + m^)/ \m e m^a) is the muonium Bohr radius.) 
The resulting mass splitting between the two mass eigenstates of the mixed M-M system 
is [34], 

AM = 2 |^| = %+^, (19) 
and the time-integrated conversion probability is 

dtF, sin 2 (AMt)e-^ = r , /(AM)2 + 4 • (20) 
The bound from the MACS experiment [33] then translates into \Y^ e + Y*^\ < 0.079. 



D. Constraints from magnetic dipole moments 



The CP conserving and CP violating parts of the diagram in Fig. 4 generate magnetic 
and electric dipole moments of the muon, respectively. Since the experimental value of the 
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magnetic dipole moment, — 2, is above the SM prediction at more than 3a, also the 
preferred value for the flavor violating Higgs couplings will be nonzero. 

The FV contribution to (g — 2)^ due to the r-Higgs loop in Fig. 4 is (neglecting terms 
suppressed by m^/m T or m T /rrih) 

q u — 2 Re(Y UT Y TU ) m u m T / , m? \ . , 

a " S V^^8^1^( 2l0g 4- 3 )' <21) 

The discrepancy between the measured value of a M and the one predicted by the Standard 
Model [30, 35], 

Aa M = a^ xp - af 1 = (2.87 ± 0.63 ± 0.49) x 10~ 9 , (22) 

could thus be explained if there are FV Higgs interactions of the size 

Re(Vy T „) ^ (2.7 ± 0.75) x 10" 3 , (23) 

(for the definition of the Yukawa couplings see Eq. (1)). This explanation of Aa M requires 
Y^ T ~ Y^ T to be a factor of a few bigger than the SM value of the diagonal Yukawa, m T /v, 
and is in tension with limits from r — > /ry 3 . It is in further tension with the LHC limit 
extracted in Sec. V of this paper. 

The measured Aa M could in principle also be explained by an enhanced flavor conserving 
coupling of the muon to the Higgs if Y^ ~ 0.15 ~ 280m At /t>. However, in this case h — > /i/i 
decays would be enhanced to a level that is already ruled out by the searches at the LHC: 
From the search for the MSSM neutral Higgs boson one obtains a bound a(gg — > h — > fifi) < 
30 x <j(gg h -> /i/i) S M or Y^ < h.hm^/v [36]. 



E. Constraints from electric dipole moments 

If the flavor violating Yukawa couplings in Fig. 4 are complex, the diagram shown there 
generates also an electric dipole moment (EDM) for the muon. The relevant term in the 



3 If the two loop contribution to t — > [i-f is suppressed, e.g. due to a modification of the top Yukawa 
coupling, which could lead to significant cancellation between the 2-loop top and W diagrams, there is 
a small region of parameter space in which flavor violating Higgs couplings could explain the (g — 2) M 
discrepancy without being ruled out by the one loop r^/J7 constraint. We will, however, see below that 
even this case is disfavored by the LHC limit derived in this paper (see Sec. VA). 
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Figure 5: Diagrams contributing to \i — > e conversion in nuclei via the flavor violating Higgs 
Yukawa couplings Y^ e and Y efl . 



effective Lagrangian is 



£edm = ~\^ (/^ Q V A*) F a p , (24) 



with the electric dipole moment given by (neglecting the terms suppressed by m tl jm T or 
m T /m h ) 



d 



'' ' 16tt 2 2m 2 



(2 log ^|-3). (25) 



The experimental constraint —10 x 10 20 e cm < < 8 x 10 20 ecm [29] translates into the 
rather weak limit -0.8 < Im^Y^) < 1.0. 

A similar diagram with electrons instead of muons on the external legs also contributes to 
the electron EDM, d e . The experimental constraint \d e \ < 0.105 x 10~ 26 e cm [29] translates 
into |Im(Y er Y re )| < l.lxl0~ 8 for a tau running in the loop, and into |Im(Y eM Y^ e )| < 9.8xl0~ 8 
for a muon running in the loop. 



F. Constraints from fi — > e conversion in nuclei 

Very stringent constraints on the FV Yukawa couplings Y^ e and Y efJi come from experi- 
mental searches for /i — > e conversion in nuclei. The relevant diagrams with one insertion of 
the FV Yukawa coupling are shown in Fig. 5. An effective scalar interaction arises already 
at tree level from the first diagram in Fig. 5, while vector and electromagnetic dipole contri- 
butions arise at one loop level. We give complete expressions for the tree level and one loop 
contributions in Appendix A3. There are also two-loop contributions, similar to the ones 
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discussed in Sec. Ill A in the context of \i — > e^y. Numerically, the two-loop contributions 
are larger than the one loop ones because they are not suppressed by the small coupling 
but only by Y tt or the gauge weak coupling, but they are still slightly less than an order of 
magnitude smaller than the tree level contribution. Here, we always assume the diagonal 
Yukawa couplings to have their SM values. With this assumption, the tree level term is very 
sensitive to the strangeness content of the nucleon. 

The bounds on the Yukawa couplings Y efM and Y^ e from /i — > e conversion in nuclei are 
listed in Table I. 

One could potentially also obtain interesting limits on \Y eT \ and \Y re \ from /i — > e conver- 
sion in nuclei, even though this requires diagrams proportional to two FV Yukawa couplings, 
because the other constraints on these couplings are weak. The combinations Y eT Y TfM , Y er Y* T , 
Y* e Y Tfl and Y* e Y* T are constrained by fi — >• e conversion through 1-loop diagrams similar to 
the ones shown in Fig. 5, but with a r running in the loop (see Eq. (A17)). In the simplest 
case, Y eT = Y Te , Y^ T = Y Tfl , with all Yukawa couplings real, the constraint is Y^Y^ < 1CT 6 . 
This is almost, but not quite, competitive with the bound following from r — > and r — > /ry 
decays, see Table I. 

G. Allowed branching ratios for lepton flavor violating Higgs decays 

In Fig. 6 we collect the above constraints on the values of \Y eT \, \Y Te \ (upper left panel), 
l^e/tl, |^te| (upper right panel) and \Y^ T \, \Y TfJ ,\ (lower panel) and relate them to the predicted 
branching ratios for h — >■ er, h — > er and h — > fir. The latter are given by 

BR(A -> «*) = r( ^V^r m • « 26 » 

where £ a , ^ = e, fi, r, £ a ^ The decay width T(h ->■ £ a £?), in turn, is 

r(fc pp) = ^(lY^l 2 + |W| 2 ) , (27) 

and the SM Higgs width is Tsm = 4.1 MeV for a 125 GeV Higgs boson [37]. In the panels 
of Fig. 6 we are assuming that at most one of non-standard decay mode of the Higgs is 
significant compared to the SM decay width. 

From Fig. 6 we see that given current bounds from r — > ^7 and r — > ej, branching 
fractions for h — > t/j, or h — > re in the neighborhood of 10% are allowed. This is well within 
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Figure 6: Constraints on the flavor violating Yukawa couplings \Y eT \, \Y Te \ (upper left panel), \Y efl \, 
\Yfj, e \ (upper right panel) and |1^ T |, \Y Tfl \ (lower panel) of a 125 GeV Higgs boson. The diagonal 
Yukawa couplings are approximated by their SM values. Thin blue dashed lines are contours of 
constant BR for h —> re, h — > fie and h — > t/jl, respectively, whereas thick blue lines are the 
LHC limits derived in Sec. VA. (These limits could be greatly improved with dedicated searches 
on existing LHC data, see Sec. VC.) Shaded regions show the constraints discussed in Sec. Ill 
as indicated in the plots. Note that g — 2 [EDM] searches (diagonal black dotted lines) are only 
sensitive to parameter combinations of the form R^Y^Ygo,) [Im(l^g ^8e*)]- We also show limits 
from a combination of g — 2 and EDM searches with marginalization over the complex phases 
of the Yukawa couplings (green shaded regions). Note that (g — 2)^ provides upper and lower 
limits (as indicated by the double-sided arrows in the lower panel) if the discrepancy between the 
measurement and the SM prediction [30, 35] is taken into account. The thin red dotted lines show 
rough naturalness limits YijYji < mirrij/v 2 (see Sec. II). 
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Figure 7: Two representative diagrams through which flavor violating Higgs Yukawa couplings can 
contribute to neutral meson mixing. 

the reach of the LHC as we shall show in Sec. V. The allowed sizes of these two decay widths 
are comparable to the sizes of decay widths into nonstandard decay channels (such as the 
invisible decay width) that are allowed by global fits [38] . If there is no significant negative 
contribution to Higgs production through gluon fusion, one has BR(h — > invisible) < 20%, 
while allowing for arbitrarily large modifications of gluon and photon couplings to the Higgs 
constrain BR(h — > invisible) < 65% [38]. These two bounds apply without change also to 
BR(h ->■ Tfi), BR(h ->■ re) and BR(h ->■ en). 

In contrast to decays involving a r lepton, the branching ratio for h — >■ e/i is extremely 
well constrained by /x — > and /i — > e conversion bounds, and is required to be below 
BR(h -> e/i) < 2 x 10~ 8 , well beyond the reach of the LHC. 

IV. HADRONIC FLAVOR VIOLATING DECAYS OF THE HIGGS 

We next consider flavor violating decays of the Higgs to quarks. We first discuss two-body 
decays to light quarks, h — > bd, bs, sd, cu, and then turn to FV three body decays mediated 
by an off-shell top, h — > i*c — > Wbc and h — > i*u — > Wbu as well as FV top decays to t — > ch 
and t — > uh. Our limits are summarized in Table II. 
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Technique 


Coupling 


Constraint 




D° oscillations [39] 


W I 2 W I 2 
[ 1 uc | > | 1 cu I 

\Y Y 1 
\ 1 uc 1 cu\ 


< 5.0 x 10~ 9 

< 7.5 x 10- 10 




B® oscillations [39] 


\Ydb\ 2 , \Ybd\ 2 
\Y db Y bd \ 


< 2.3 x 10~ 8 

< 3.3 x 10- 9 




oscillations [39] 


\Ysb\ i \Y~bs\ 

\Y sb Y hs \ 


< 1.8 x 10~ 6 

< 2.5 x 10~ 7 






Re(Yl), Re(Y s 2) 


[-5.9... 5.6] x 10- 


-10 


K° oscillations [39] 


im(^). Im (n 2 d ) 

MY d * s Y sd ) 


[-2.9... 1.6] x 10" 
[-5.6... 5.6] x 10" 


-12 
-11 




MY d * s Y sd ) 


[-1.4... 2. 8] x 10- 


-13 


single-top production [40] 


Vl^el + 
vWl + \Yut? 


< 0.54 

< 0.23 




t->hj [41] 


^\Yl\ + |Y*| 2 
y/X*l\ + \Y ut \ 2 


< 0.34 

< 0.34 






\Y u tY c t\, \Y tu Y tc \ 


< 7.6 x 10~ 3 




D° oscillations [39] 


\Yt u Y c t\, \Y ut Y tc \ 
lYutYtuYctYtd 1 / 2 


< 2.2 x 10~ 3 

< 0.9 x 10-3 




neutron EDM [29] 


Im(Y ut Y tu ) 


< 4.4 x 10~ 8 





Table II: Constraints on flavor violating Higgs couplings to quarks. We have assumed a Higgs mass 
vrih = 125 GeV, and we have taken the diagonal Yukawa couplings at their SM values. 

A. Flavor violating Higgs decays into light quarks 

Flavor violating Higgs couplings to quarks can generate flavor changing neutral currents 
(FCNCs) at tree level, see Fig. 7 (a), and are thus well constrained by the measured B djS — 
B djS , K° — K° and D° — D° mixing rates. Integrating out the Higgs generates an effective 
weak Hamiltonian, which for B d — B d mixing is 

H cS = C* b (b R d L ) 2 + Cf{b L d R f + Cf(b L d R )(b R d L ) . (28) 
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Here we use the same notation for the Wilson coefficients as in [39] and display only nonzero 
contributions, which are 



(V*\2 (V2\2 V V* 

n db _ Vdb) f<db _ \ I bd) n db _ I bdidb / on \ 

2 ~~^P 2 ~~W 4 ~~^T' { } 



The results for B s — B s , K° — K° and D° — D° mixing are obtained in the same way with the 
obvious quark flavor replacements. We can now translate the bounds on the above Wilson 
coefficients obtained in [39] into constraints on the combinations of flavor violating Higgs 
couplings as summarized in Table II. We see that all Yukawa couplings involving only u, 
d, s, c, or b quarks have to be tiny. The weakest constraints are those in the b-s sector, 
where flavor violating Yukawa couplings < 10~ 3 are still allowed. This would correspond to 
BR(/i — > bs) ~ 2 x 1CT 3 , which is still far too small to be observed at the LHC because of 
the large QCD backgrounds. 

B. Higgs decays through off-shell top and top decays to Higgs 

Among the flavor violating Higgs couplings to quarks, the most promising place for new 
physics to hide are processes involving top quarks, such as the 3-body decay h — > (t* — > 
Wb)q. Here, q denotes either a charm quark or an up quark. The corresponding FV Yukawa 
couplings contribute at one loop to D — D mixing through diagrams of the form of Fig. 7 (b). 
The corresponding Wilson coefficients in the effective Hamiltonian (28) are 

Cl "4 16^ ml ( YctY ^ ' Cl ~lW m\ ' (3 ° } 

° 2 - _ 416vr2 ml [YtcKt) ' ° 2 "~416tt2 ml ^ ctYtu) ' (31) 

cr = -~/ 1 ^(Y ct Y tu )(Y t : Y : t ), C r = - i ^^)(w)(n^), (32) 

where 

qH / \ _ a: 2 -1 -2a: log x „ _ 2a;[2 - 2x + (1 + x) logs] 

bi[x) ~ 2(x-i)3 ' 2( j ~ (x~^rr ' ( j 

and x t H = ml/ml. Note that now also the operators Qi^, Q\ c (in the notation of [39]) have 
non-zero Wilson coefficients. By requiring that each individual operator is consistent with 
its D — D mixing constraint, we derive the limits shown in the last part of Table II. The 
constraints are much weaker than those on FV Higgs couplings involving only light quarks. 
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Figure 8: Predictions for various flavor changing neutral current (FCNC) processes mediated by the 
flavor violating Yukawa couplings Yet, Yt c or Y u t, Ytu of a 125 GeV Higgs boson. Where appropriate, 
we have approximated the diagonal Yukawa couplings by their Standard Models values. Blue 
dashed contours indicate the branching ratio for h — > t*q, red solid contours the one for t —> hq 
(where q denotes a charm or up quark). The red dotted line is a recent limit on t — > he (or hu) 
from an LHC multi lepton search [41]. 

Strong constraints on Y qt and Y tq are also obtained from the non-observation of anomalous 
single top production. The flavor violating chromomagnetic operators 

9 \ a 

^single top ^ (K tqg L P L + K, tqgR P R )—q G°^ v , (34) 

771/j L 

are generated trough loop diagrams similar to Fig. 1, but with leptons replaced by quarks 
and the photon replaced by a gluon. Here g s is the strong coupling constant, A a are the Gell- 
Mann matrices, G a is the gluon field strength tensor, and Ktqg,L, ^t qg ,R ar e dimensionless 
effective coupling constants which depend on Y qt and Y tq according to 



K — sb^(- 4+31 °^)- < 35 > 



The analogous expression for Kt qg ,R is obtained by replacing Y t * — > Y qt and Y tt — > Y t *. 
Limits on Ktqg,L, K tqg,R have been derived by the CDF and D0 collaborations [40, 42] 
and most recently by ATLAS [43]. In the notation of [43], we have |^ 9 /|/A = 
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\/\ K tqg,L\ 2 + \K-t qg ,R\ 2 / (V2rn h ). We obtain the constraints 



^/|y t 2| + |y rf |2<o.54, ^| + |y ut |2<o.23, (36) 

We now translate these bounds into constraints on the h — > (i* — > decay width, which 

is given by (setting m^g = 0) 

d 2 T(h^i*q)_ 3g 2 \V tb \ 2 1 
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(2m 2 w - my (m 2 ^ 2 - m 2 23 \Y tq \ 2 ) 



dm\ 2 dm% 3 6A(2Tv) 3 m 2 v m 3 h (m| 3 - m?) 2 

+ K - m 2 3 ) ("4 - "4) (2m^|y t? | 2 + m 2 \Y qt \ 2 )^ , (37) 

where V th ~ 1 is a CKM matrix element. The branching ratio for /i — >■ t*c can be as large 
as a few per cent, and the one for h — > t*u a few per mille as shown in Fig. 8. 

If the decay h — > (t* — > Wb)c is non- negligible, so is the related non-standard top quark 
decay mode t — > he, the rate for which is given by (neglecting the charm mass) 

rNk) Jtf + ltfW T D' (38) 

v ; 32tt mf y J 

Branching ratios for t — > he of several tens of per cent are perfectly viable and can be 
searched for, e.g. in the multi-lepton or t — > bbc channels. In fact, the strongest hint on 
Higgs couplings to tc are already coming from a CMS multi-lepton search which was recast 
in [41] to search for t — > he, giving a bound of 2.7% on the branching fraction of a top into 



a Higgs and a charm or up quark. This yields a limit of a/I^I 2 + \Y it \ 2 < 0.34 for % = u or c 
(see Fig. 8). 

We have also calculated the branching ratios for the loop-induced processes t — > q^y, 
t — > qg and t — > qZ (q = u,c), which are in principle sensitive to \Y qt \ and \Y tq \, but have 
found that even for \Y qt \, \Y tq \ ~ 0(1) the current experimental bounds are satisfied [44]. 

In the above we have assumed that the weak phases of Y ut and Y ut are negligibly small. 
Otherwise an unacceptably large contribution to neutron the EDM is generated at 1-loop 
level with top and Higgs running in the loop. Eq. (25) with m T — > m t , and Y^ T Y T ^ — > Y ut Y tu 
replacements gives the w-quark EDM d u = ed u , from which one can calculate the neutron 
EDM d n [45]. Using the 90%CL experimental bound d n < 0.29 x 10~ 25 ecm [29] together 
with the theoretical prediction, Eq. (3.62) of [45] one obtains |Imy ut y tu | < 4.4 x 10 -8 , 
which is much more stringent than the bounds on the absolute values of the same FV 
Yukawa couplings. In contrast, the bounds from charm running in the loop, |ImY Mc y cu | < 
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1.6 x 1CT 7 , and from d-quaik EDMs generated by the 6-quark and s-quark running in the 
loop, | Im Y^Y^ | < 6.4 x 1CT 8 and |ImF rfs F sd | < 1.2 x 10~ 6 , respectively, are less stringent 
than the bounds from meson mixing, Table II. 

V. SEARCHING FOR FLAVOR VIOLATING HIGGS DECAYS AT THE LHC 

We next discuss possible search strategies for flavor violating Higgs decays at the LHC, 
focusing on the h — > t/i and h — >■ re decays. As shown in Fig. 6, these are among the least 
constrained of the couplings discussed in this paper, with a potential to modify the Higgs 
branching fractions significantly. They are sensitive to new particles with flavor violating 
couplings or to a secondary mechanisms of electroweak symmetry breaking such as addi- 
tional Higgs doublets, and are thus good probes of new physics. Furthermore, they are also 
interesting final states as far as the potential for searches at the LHC is concerned. 

The decay h — > t/i is quite similar to the standard model h — > tt decay with one of 
the tau leptons decaying to a muon. This implies that existing SM Higgs searches, with 
only small or no modifications at all, can already be used to place bounds on the flavor 
violating decay. We thus first extract limits on h — >■ t/i and h — >■ re decays from an existing 
h — > tt search in ATLAS. We then discuss how modifications to the tt search can lead to 
significantly improved sensitivity to flavor violating Higgs decays. 

A. Extracting a bound on Higgs decays to r/z and re 

We use the existing ATLAS search for h — >■ tt in the fully leptonic channel [46] to place 
bounds on the h — > t/i and h — >■ re branching fractions. The reason we use fully leptonic 
events is that we can simulate the detector response to them more accurately than for events 
involving hadronic taus. It should, however, be noted that in the SM h — >■ rr search, semi- 
hadronic events are about as sensitive as fully leptonic ones [46]. The analysis in [46] uses 
the collinear approximation to reconstruct the rr invariant mass, i.e. it is assumed that the 
neutrino and the charged lepton emitted in tau decay are collinear. This approach is less 
optimized for h — > rr than the maximum-likelihood method employed by CMS [47], but it 
is more model independent so that a substantial fraction of h — > t/i or h — > re decays would 
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pass the cuts 4 . For simplicity we only use the ATLAS cuts optimized for Higgs production 
in vector boson fusion ( VBF) since this channels provides the best sensitivity [46] . 

To derive limits we have generated 50,000 pp — > 2j + (h — > t/j) Monte Carlo events using 
MadGraph 5 vl.4.6 [48] for parton level event generation, Pythia 6.4 for parton showering 
and hadronization, and PGS [49] as a fast detector simulation. Combining the ATLAS 
lepton triggers and off-line cuts from [46], we select opposite sign dilepton events satisfying 
any of the following requirements: a muon pair with p T > 15 GeV for the leading muon and 
p T > 10 GeV for the subleading one, an electron pair with both p T > 15 GeV, or an electron 
and a muon with p T above 15 and 10 GeV, respectively. Electrons (muons) are accepted 
only if their pseudorapidity is \r)\ < 2.47 (2.5). We require the invariant mass of the lepton 
pair to be 30 GeV < mg < 100 GeV for e/x pairs, or 30 GeV < ran < 75 GeV for same flavor 
pairs. The missing px is required to be above 20 (40) GeV for e/i events (ee or events). 
The azimuthal separation between the two leptons is required to be 0.5 < A(pn < 2.5. 

Additional cuts are placed with the goal of enriching the event sample in VBF events: at 
least two jets with p T above 40 GeV for the leading jet 25 GeV for the subleading jet are 
required, with the rapidity difference between the two leading jets above \Ar]\ > 3 and the 
rapidity difference rrijj > 350 GeV. We veto events with an additional jet with pj- > 25 GeV 
and \r)\ < 2.4 in the pseudorapidity region between the two leading jets. 

The reconstructed invariant mass is calculated using the collinear approximation in which 
all invisible particles are assumed to be collinear with either of the two leptons. The fractions 
of the parent r's momenta carried by the charged leptons are denoted by x\ and x 2 . To be 
able to compare with ATLAS data from the h — >■ tt search, we compute x\ and x 2 assuming 
two neutrinos in the final state, even though h — > t/j, yields only one. x\ and x 2 are then 
obtained as the solutions of the transverse momentum equation p m i SS ,T = (1 — ^i)Pi,t + 
(1 — x 2 )p2,t, where Pi,2,t are the transverse momenta of the charged leptons. Following [46], 
we require 0.1 < x 1>2 < 1, which removes less than a per cent of h — >■ tt events, but nearly 
60% of our h — >■ Tfi events. Thus, relaxing this cut would enhance the sensitivity to h — > t/j 
decays so long as it does not introduce large backgrounds. Nonetheless, we are still able to 

4 In fact, it may be interesting to apply the collinear approximation more often in resonance searches. A 
search for a collinear mass resonance can be sensitive to any particle which decays to boosted objects 
whose further decay may introduce missing energy. 
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Figure 9: Background rates and h — > Tfi, h —> tt signal rates in the ATLAS search for fully 
leptonic h — > tt decays, optimized for Higgs production in vector boson fusion. The backgrounds 
expected by ATLAS [46] are shown in yellow, with grey bands for the systematic uncertainty. Our 



estimates for the t\i signal at + Yj* T = m T /v (red) and the SM h —> tt signal (black), which 

we include for reference, are scaled by a factor 5 for illustrative purposes only. 

use the current search for h tt to produce an interesting bound on BR(/i — > t(j). 

In Fig. 9 we show the background distribution for the collinear mass along with the 
expected shape of a LFV h — >■ r/x signals (scaled by a factor five for illustrative purposes 
only), and we compare to the observed data. The background expectation is taken from [46]. 
The backgrounds and the data in Fig. 9 includes events for all three combinations of lepton 
flavor (even though our r\x signal does not induce ee events) because only this information is 
available from ATLAS. For validation purposes, we have also simulated SM h — » rr events, 
and comparing the rate and shape to Ref. [46] we find agreement to within 20%. 

The t/j, signal is predominantly concentrated in the 120-160 GeV bin, so that the expected 
and observed limits on the flavor violating Yukawa couplings can be derived from a simple 
single-bin analysis. If we denote the number of expected background events by B = 4.7, the 
number of expected signal events for a given set of Yukawa couplings by S, and the number 
of observed events by O = 2, the expected (observed) one-sided 95% CL. frequentist limit 
on S is defined by the requirement that the probability to observe < B (< O) events is 
5%. The relevant probability distribution of the data here is a Poisson distribution with 
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95% C.L. limit 


BR(h -»• Tfi) 




BR(h -»• re) 




expected 


28% 


0.018 


27% 


0.017 


observed 


13% 


0.011 


13% 


0.011 



Table III: Expected and observed 95% C.L. limits on the h — > tu and h — > re branching ratios, 
as well as limits on the corresponding Yukawa couplings. The limits are derived by assuming the 
SM Higgs production rates and recasting the search for SM h — > tt — > 21 + 2v decays in the VBF 
channel from [46]. 

mean B + S. We can also include the systematic uncertainty in the 120-160 GeV bin, 
which is A sys ~ ±0.99, in a conservative way by instead using a Poisson distribution with 
mean B + S — A sys . Assuming the Higgs is produced with the Standard Model rates, this 
procedure leads to the bound on BR(/i — > t/j) and the analogous bound on BR(h — > re) 
shown in Table III (see also Figure 6). 

B. Comparison of h — >■ r/x to h — ^ rr 

We now discuss the experimental differences and similarities between h — >■ rr and h — > r[i 
decays to determine an optimized search strategy for the latter. We focus here on h — > 
ThadT/j,, where r M denotes a r that decays into a muon and two neutrinos and rh a( j denotes 
a r decaying hadronically. This channel is actively searched for, both at ATLAS [46] and 
at CMS [23], and is the most sensitive channel in the CMS h — ?> rr search. (In ATLAS, 
fully leptonic r events provide similar sensitivity to semi-hadronic ones.) It will also be the 
channel that we will devise a dedicated search for in the next subsection. 

There are a few notable differences between the h — > ri ia dr M and h — > rh a d A* decay channels: 

• Branching Ratios. The branching fraction for h — > r had r M is 2 x BR(/i — > rr) x BR(r — > 
had) x BR(r — > //), whereas for h — > r had /i it is simply BR(h — > t/i) x BR(r — > had). 
For Y Til ~ Y TT the signal for h — > r ha dAt is thus a factor of ~ 2.9 larger. 

• Lepton Flavor. The flavor violating decays can lead to different rates for muons and 
electrons in the final state, whereas rr decays lead to equal /i and e rates. Thus, if the 
various lepton flavor combinations were studied separately in the h — > rr analyses, 
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stronger bounds on flavor violating decays could be inferred. 

• Kinematics and Efficiencies. In h — > Th ad T M decays the muon carries an average energy 
~ rrih/6, while for h — > Th ad /i it carries ~ m/j/2. Furthermore, in h — > Th ad H events the 
missing energy is roughly aligned with the hadronic r. As a result the two channels 
will have different efficiencies given the same cuts. For example, in the VBF analysis 
described below (mimicking [23]) the efficiency for h — > r had r M is a factor of ~ 2.3 
lower than for h — > Th ad/ u events, mostly because many of the muons in the h — > t^&t^ 
sample fall below the pr < 17 GeV cut. 

• Mass reconstruction. The LHC collaborations use highly optimized procedures for 
reconstructing the Th ad T M invariant mass. ATLAS uses the Missing Mass Calculator 
(MMC) from [47], while CMS uses an in-house maximum likelihood analysis [23]. 
These procedures use p m i SS ,T and the 3-momenta of the muon and the r jet as input 
and estimate the neutrino momenta by assuming typical r decay kinematics. For 
h — )• tt events, the MMC procedure returns an invariant mass with high efficiency 
(~ 97%) and gives a Higgs mass resolution of ~ 20%. If the event is not from h — >■ tt 
but instead from h — > Th a d /i, then i) the efficiency will be significantly lower since the 
kinematics can be completely inconsistent with a tt event, and ii) the reconstructed 
Higgs mass will be significantly higher as the MMC will assume that the hard muon 
is accompanied by two roughly collinear and hard neutrinos. This illustrates that a 
mass reconstruction procedure designed for the specific final state under consideration 
is mandatory to obtain the best possible sensitivity. 

• Backgrounds. The backgrounds for h — > r had r M and h — > r had fj, events are similar, but 
because of the different invariant mass reconstruction techniques, the reconstructed 
background spectra will typically be harder for a h — > Th a dT M analysis, which assumes 
four neutrinos in the final state, than for a h — > Th a d/i analysis which assumes only one. 
This implies, for instance, that the peak from the Z — >■ tt background will appear at 
a T had r M invariant mass around 90 GeV in a search for h — > t^^t^, but well below (and 
thus further away from the signal peak) in a dedicated h — > Th ad /i analysis. 

These considerations show that the LHC is potentially more sensitive to flavor violating 
h — > T had /x decays than to the SM h — >■ tt channel. We now discuss a possible strategy for 
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a tailored h — > Th a d/i analysis. 

C. A dedicated h — > t\i analysis 

We now investigate the potential of a dedicated h — > r ha dAt analysis which follows closely 
the CMS search for h — > twit^ [23]. The most important difference to that analysis will be 
a different algorithm for reconstructing the t\i invariant mass. In particular, since the Th a d/^ 
final state contains only one neutrino (from the hadronic r), this mass reconstruction can 
always be done exactly (i.e. the neutrino momentum can be determined) up to a two-fold 
ambiguity. 

An important background for h — >■ t/i is Z + jets, where either the Z decays into t + t~ 
and one of the r's decays further into a muon, or the Z decays into fi + {i~ and one of the 
jets fakes a r. Another important background is W + jets, followed by W — > jiu^ and a jet 
faking a r. We neglect the small ti background, where a final state r can come from a W 
decay or be faked by a jet, and a muon can originate from a W decay or from a leptonic r 
decay. We also do not consider backgrounds from QCD multijet production because making 
reliable predictions for these events requires full detector simulations. Based on the CMS 
h — > r M r had search [23] we expect them to be about as large as the W + jets background, in 
the invariant mass region around 125 GeV. 

To simulate the parton-level signal and background events, we use MadGraph 5 vl.4.6 
[48], with an extended version of the Higgs Effective Theory model to include flavor- 
violating Higgs interactions. We use Pythia 6.4 for parton showering and hadronization 
and Delphes 2.0.2 [50] as a fast detector simulation. We have compared the r detection 
efficiency as well as the fake rate from QCD jets in Delphes 2.0.2 to the corresponding per- 
formance indicators of several CMS r tagging algorithm. With a tagging efficiency of ~ 0.2 
and a fake rate between 0.2% at low pr and 1% at pr > 60 GeV, Delphes somewhat under- 
estimates the performance of the CMS r tagging algorithms [51]. (We have also studied r 
tagging in PGS [49], but found it to be even farther away from what CMS can achieve.) To 
compensate for the imperfections in our treatment of r-tagging, and for other inaccuracies in 
our simulations, we normalize our background distributions to the expected event numbers 
from [23], Table 2. We also normalize the h — > r/i signal using the same scaling factor as 
for the SM h — > events. 
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In the analysis we require exactly one muon with px > 17 GeV and \r]\ < 2.1 in the final 
state and exactly one jet tagged as a hadronic r decay with p T > 20 GeV and \r)\ < 2.3. 
The muon and the r are required to have opposite charge. In [23], it was found that 
the best signal-to-background ratio is achieved in the events where the Higgs boson was 
produced through vector boson fusion (VBF), and we confirm this in our own simulations. 
To enrich the data sample in VBF events, we consider only events with a pair of jets ji, j2 
satisfying \Ar/\ > 4.0, 771772 < 0, rrijj > 400 GeV and no other jets with p T > 30 GeV in 
the pseudorapidity region between the ji and j 2 . Here, At] = 771—772 is the pseudorapidity 
difference between the two jets and rrijj is the invariant mass of the jet pair. Non-r jets 
are included in the analysis so long as their pt is above 30 GeV and their pseudorapidity 
is 1 7/| < 4.7. In the CMS analysis [23], the transverse mass of the muon and the missing 
energy is restricted to be below 40 GeV in order to suppress the W + jets background. This 
works because in W + jets events, the muon and the neutrino from W — > /iv^ tend to be 
more back-to-back than in h — > r M r had , where both r's contribute to the missing energy. In 
h — > T had /i, however, the muon and the missing energy also tend to be back-to-back, so that 
the mxifi, Pmiss.T) also removes a large fraction of signal events. 

In light of this we show in Fig. 10 the expected signal and background rates for h — > t/j, 
as a function of the /i-r invariant mass m Tll both with and without the transverse mass 
cut. In computing m T ^ for each event, we have used energy and momentum conservation 
to compute the ^-component of the neutrino momentum p vz . There are two solutions to 
these equations, and we arbitrarily pick the smaller of the two. (We have checked that 
choosing the larger value for p U:Z yields a very similar plot. This is related to the fact that 
m T <C rrih, so that the r's decay products are almost collinear.) As shown in the right 
panel of Fig. 10, dropping the transverse mass cut increases the W plus jets background, 
but has the benefit of retaining more signal. The transverse mass distributions for signal 
and background is shown in Fig. 11. Assuming the W plus jets and QCD backgrounds can 
be controlled reasonably, relaxing this cut may be worthwhile. 

In summary, Fig. 10 shows that for flavor violating Yukawa couplings well allowed by low 
energy precision measurements, a spectacular signal can be expected in a dedicated search 
at the LHC. Such a search would cut deeply into the allowed parameter space of the flavor 
violation Higgs to r/i couplings. 
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Figure 10: Signal and background rates for h — > r/x events in a CMS-like search (see text) as 
a function of the reconstructed \x-t invariant mass m TfJi for a vector boson fusion-enriched event 
sample. In the left panel the transverse mass cut m^di, p m iss,T) < 40 GeV is included, while in 
the right panel it is omitted. The QCD multijet background and the small tt background, are not 
included. The value chosen for <Jy^ + Yff T is well within the region allowed by other searches for 
flavor violation in the fi-T system (see Sec. III). 
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Figure 11: The muon transverse mass distribution for the backgrounds and for the ThadA* signal. 
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VI. CONCLUSIONS 



The LHC experiments have recently discovered a Higgs-like resonance with a mass around 
125 GeV. In this paper we have examined the constraints on potential flavor violating 
couplings of this resonance, assuming it is indeed a scalar boson. In deriving the constraints 
we have assumed that that flavor changing neutral currents are dominated by the Higgs 
contributions, which may be thought of as a "simplified model" approach to flavor violation 
in light of the Higgs discovery. (In a complete model, cancellations between Higgs-induced 
flavor violation and flavor violation induced by other new physics is possible, but we do not 
pursue this possibility here.) 

We have derived new constraints on the flavor violating Yukawa couplings Yy using results 
from rare decay searches, magnetic and electric dipole moment measurement, and the LHC. 
All constraints are summarized in Tables I and II and in Figs. 6 and 8. We have compared 
the bounds to the loose naturalness criterion that the off-diagonal Yukawa couplings are 
not much bigger then the geometric mean of diagonal terms, < y / Y^Yjj, and we have 
discussed to what extent the LHC can probe flavor violating decays of the form h — > fifj 
(or in the case of flavor violating top-Higgs couplings the decay t — > hfi). 

We draw the following conclusions: 

Natural Flavor Violation. The existing constraints involving only the first two gener- 
ations of fermions, quarks or leptons, are strong enough that natural FV is already being 
probed by meson oscillations and /i — > e conversion. This conclusion also holds for FV 
couplings to b quarks. In contrast, the FV couplings involving r leptons or top quarks 
are allowed to have natural size (unless there is a large hierarchy between and Yji). 
This means that they are potentially observable, either at the LHC or in future low energy 
experiments. 

Opportunities for the LHC. The LHC has an opportunity to probe a large part of the 
allowed parameter space for h — > r/i and h — > re couplings. An LHC search in these 
channels would be very similar to the existing searches for h — > rr, and recasting the latter 
already gives the best bounds on the flavor violating Yukawa couplings Y re , Y er , Y TfJi , and 
Y^ T already now. A dedicated LHC search could improve the limits significantly. The reason 
why Higgs decays are very constraining is that the SM width of a 125 GeV Higgs boson is 
very small, Y H ~ 4 MeV [37], so that the flavor violating couplings of the Higgs can have a 
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significant effect. Another illustration of the LHC's discovery potential for flavor violating 
Higgs couplings is that even the global fits of potential deviations in the dominant SM Higgs 
decay modes, h — > WW, ZZ, bb, tt, 77, already give meaningful bounds on the FV Higgs 
decays. The results from these global fits are usually presented as bounds on the invisible 
decay width of the Higgs, but these bounds applies equally well to the sum of all the modes 
that have not been included in the fits. The constraint BR(h — y invisible) < 0(70%) (at 
95 % CL, with modest theory assumptions [38, 52-54]) is comparable to the constraints on 
BR(h — y Tfi) and BR(h — y re) from precision searches of FCNCs in the lepton sector. 

Finally, flavor violating Higgs couplings involving the top quark are poorly constrained, 
and in fact, the LHC already provides the strongest limits on such couplings, see Fig. 8 and 
Section IV. 
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Appendix A: Further details on leptonic FCNCs 

In this appendix we collect detailed expressions for the FCNC processes r — y /ry and 
fj, — y e conversion in nuclei. 

1. One loop expressions for r — y ^7, r — y ej, ji — y e'y 

The t — y /i^f effective Lagrangian is given in Eq. (11). The Wilson coefficients cl,r are 
given by 

c L = F(m T , m T , ra M , m h , 0, Y) + F(m T , m M , m M , m h , 0, Y) , (Al) 
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c R = F(m T , m T , m M , m h , 0, Y r ) + F(m T , m M , m M , m h , 0, Y f ) , (A2) 
with the loop functions 

1 

dx dy dz 5(1 — x — y — z) 

xzm 3 Y jf Y* f + yzmjYfjYfi + (x + V^fY^Y* (A3) 
zm\ — xz m 2 — yz mf + (x + y)m 2 — xyq 2 

Here m M , m T and are the muon, tau and Higgs masses, respectively, q is the 4- momentum 

of the photon, and Y is the Yukawa coupling matrix. Note that cl and cr differ only by 

the replacement Y^- -H- Y^*. The first terms in Eqs. (Al) and (A2) arise from the first 

diagram in fig. 1 (with a r propagating in the loop), whereas the second terms arise from 

the second diagram (with a // in the loop). Expanding in powers of m M /m r and m r /m/ l 

and keeping only the leading terms (so that only the first terms in (Al), (A2)), the above 

expressions simplify to (13) for real diagonal Yukawa couplings. The simplified expressions 

for r — > e7 and /i — > cy (with a muon running in the loop) are obtained from (13) with 

trivial modifications, while the simplified expression for /i — > ey with a r running in the loop 

is given in Eq. (16). 



F(m h m s ,m^m h ,q 2 ,Y) = 




2. Two loop expressions for r — > fij, r — > e^f and \i — > 

At two loops there are numerically important diagrams with top or W running in the 
loop, attached to the Higgs. Here we translate the results of [55] into our notation and adapt 
them to the case of r — > /ry. The diagrams with top and photon in the loops (see Fig. 12 
top left) contributes as 

A$ = -QkQI—Y^ [Re{Y u )f{z th ) - ilm(Y tt )g(z th )] , (A4) 

while the IT^-photon 2-loop contribution is 



Sf(z Wh ) + 5g(z Wh ) + \g{z Wh ) + \h{z Wh ) + f{Zwh ] 9{Zwh) ' . (A5) 

IZwh 



Here we have already added the contributions from the would-be Goldstone bosons that 
get eaten by the W. The contributions to Ac R are obtained from the above by replacing 
Y*^ — > Y^ and Y a — > Y^.The loop functions are 

,/ n 1 f 1 i l-2x(l-x) , x(l -x) /A . 

f(z) = -z / dx \ >- log ^ >- , A6 

2 J x(l — x) — z z 
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Figure 12: The two loop diagrams contributing to r — >■ fij. 



h{z) = z 2 



-z [ dx— 

2 J x(l-x)-z 

i 



log 



g /g(g) 



xfl — x) 



2 J Q z — x(l — x) 



z — x(l — X 



log 



x(l — x) 



the arguments are z t h = m 2 /m H , zy/h = m w/ m Hi while the prefactor is 



K 



a g v 



a v 



167T mly m T 2y/2n m T ' 

The contributions from the 2-loop diagrams with an internal Z are smaller as 
suppressed by 1 — As w ~ 0.08. They are 

{l-4s 2 w )(l-AQ t s 2 w ) v 



(A7) 
(A8) 



(A9) 
they are 



Acf 



AcY z 



-QnQ t 



1«„2 2 



Y*.x 



m t 



K 



x [Re(Y tt )f(z th , z tz ) - ilm(Y tt )g(z th , z tz )] , 
i ■ Y *A |( 5 - tw)fi z thi z wz ) + 1(7 - 3t 2 v )g(z t h, z wz ) 



1 - 4s 2 w 



+ \g{zth) + |M^«») + 4^( x - f w) [f( z th, zwz) - g(z t h, z W z)] } , 

with s w = sin 6^, c w = cos9 w , tw = tan 6^, £t 2 = m 2 /m 2 z , z W z = ^wl m z an d 



(A10) 

(All) 
the loop 



32 



functions 

f {x ,y) = >>m + ?m, g (x ,y ) = y^l + ^. (A12) 

y — x x — y y—x x — y 

The Ac R are obtained by replacing Y* e — > Y eii and Y a — > Y£ in the above expressions. In 
addition there are also contributions called "set C" in [55] diagrams. An example for one of 
these diagrams is the last diagram in Fig. 12. Since the expressions for these diagrams are 
long we do not write them out explicitly. The "set C" contribution to Acl is obtained from 
[55] by multiplying their Eq. (20) by — k and replacing J2 a COSi Pa,A^ — > Y*^. 

The r — > expressions are obtained by replacing Y Tfl — > Y Te and Y^ T — > Y er in the above 
expressions, while for \i — > the replacements are Y T ^ — > Y^ e , Y^ T — > Y eii and m T — > m M . 



3. Details on fi — > e conversion bounds 

The most general effective Lagrangian for jj, — )■ e conversion in nuclei is [56] 
C = c L^rn^{e(T a(3 P L fi)F a p - ^ ^ [^L(e-PflA*)(99) + 9 9 Lp(eP R n)(q^ 5 q) 



q 

+ L o R. 



(A13) 



The above operators are generated by integrating out the Higgs in the diagrams shown in 
Fig. 5. The magnetic penguin Wilson coefficients Cl and cr are (in analogy to Eqs. (Al) 
and (A2)) 

c L = F(m M , m M , m e , m h , -ml, Y) + F(m jU , m e , m e , m h , -mj, Y) , (A14) 
c R = F(m^ m M , m e , m h , -mj, F f ) + F{m^ m e , m e , m h , -mj, F f ) , (A15) 

with the function F(mi,mf,mj,mh,q 2 ,Y) given in Eq. (A3). Here we have set q 2 ~ — m 2 
for the 4-momentum transfer to the nucleus, which is justified in the limit of an infinitely 
heavy nucleus. The above contributions are obtained from the last two diagrams in Fig. 5 
with either a muon or an electron running in the loop. These diagrams also generate the 
quark vector current operators with (in general gf^ = 2g^ v + gf v and g^y = ^giv + 9lv) 

9lv = 9 q LV /Qg = -^[ G K m f" m e> m fo9 2 - y ) - "V, m c , m h , 0, y) 

Z7T q z 
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+ G(m M , m e , m e , m h , q 2 , Y) - G{m^ m e , m e , m h , 0, Y)] . (A16) 

Here, Q q is the charge of quark q, and g RV , g RV is given by Eq. (A16) with the replacement 
Y ->■ F f . The loop function in (A16) is 

G(mi,m f , mj,m h , q 2 , Y) = J dx J dy Y jf Y* f log A - ^mimjZ 2 Y^Y n 

- \ ( m f m J zY fj Y *f + m f m i zY jf Y fi + W x v + m )] Y if Y ?f) , 

(A17) 

where we have defined A = zm 2 h — xz m 2 — yz m 2 + (x + — a;yg 2 and z = 1 — x — y. Note 
that we subtract the value of the one-loop vertex correction at q 2 = 0, which gets absorbed 
into the wave function and mass renormalizations. 

The coefficients cl,r also receive two loop contributions with either a top or a W running 
in the loop. If we evaluate these contributions at q 2 = 0, so that the expressions are exactly 
the same as in the previous subsection on /j, — > decay, we find that both one-loop and 
two-loop contributions are subdominant (with the 2-loop contribution the larger of the two). 
They are roughly an order of magnitude smaller then the scalar contributions, so that the 
above approximations suffice for our purposes. 

The scalar Wilson coefficients are generated from the first diagram in Fig. 5 

9ls = --V; e Re(>g , gl s = -^ Y e,Re{ Y ^ > ( A18 ) 

while the remaining Wilson coefficients are zero, g q LP = g q RP = g q LA = g q RA = g q LT = g q RT = 0. 

When computing the n — > e conversion rate in the nuclei care must be taken to account 
for the nuclear matrix elements (N\qq\N), (Nlq^qlN) and (NlF^lN) and for the overlap 
of the initial muon wave function and the final state electron wave function. We follow [56] 
and obtain 



e conversion) =| - j^c L D + g^S® + g^S™ + g^V^ 
+ ~ ^- 2 crD + ~g<&SW + g^ + 



(A19) 



16tt 2 

The electromagnetic penguin and vector contributions are given in Eqs. (A14), (A15), (A16), 
while the scalar coefficients are 

9ls,rs = ^2 9ls,rs Zt/ (9 ' p) ' 9ls,rs = ^2 9ls,rs ' (A20) 

q y q y 
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and the sum runs over all quark flavors, q — u,d, s, c, b, t. Note that vector couplings to 
neutrons are absent due to the neutron's vanishing electric charge. 

The nucleon matrix elements f^'^ = (p\m g qq\p) / m p are calculated according to [57], but 
using an updated value for the nucleon sigma term E^jv = 55 MeV [58] (If the value of the 
nucleon sigma term is even smaller, as indicated by recent unquenched lattice results, our 
bounds would become weaker). The nucleon matrix elements are numerically 



f(u,p) = f(d,n) = q Q24 ? j{d,p) = j{u,n) = Q^gg ^ j{s,p) = j{s,n) = ^ 



with the same values for neutron. In the above expressions, m q denotes a quark mass, m p 
is the proton mass, and m n is the neutron mass. The coefficients D, S^ p \ and 

are overlap integrals of the muon, electron and nuclear wave function. They are tabulated 
for various target materials in [56]. The best limits are obtained from bounds on ji — )■ e 
conversion on gold, T(fi — > e^u/Tcapturc a u < 7 x 1CT 13 (90% CL) [59], for which in units 
of mj 2 the overlap integrals are D = 0.189, = 0.0614, = 0.0974, = 0.0918, 
y(n) _ Q_x4g ; us ing the same distributions for neutrons and protons in the nucleus. For the 
SM capture rate, we use a value r captU re Au = 13.07 x 10 6 s -1 in the calculation [56]. 
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